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^— n . Abstract 

This paper is addressed to an inverse stochastic hyperbolic problem with three unknowns, i.e., a 
random force intensity, an initial displacement and an initial velocity. The global uniqueness for 
this inverse problem is proved by means of a new global Carleman estimate for the stochastic 
hyperbolic equation. It is found that both the formulation of stochastic inverse problems and 
the tools to solve them differ considerably from their deterministic counterpart. 
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> ■ 1 Introduction 

o: 

CO 



X 



Let T > 0, and let G C M. n (n £ N) be a given bounded domain with a C 2 boundary T. Put 

A ,„ _ N „ , _ A 



CO I Q = (0, T) x G and E = (0, T) x T. Fix a complete filtered probability space (fi, J 7 , {.Ft}t>o, -P), on 

f — ■ which a one dimensional standard Brownian motion {B(t)}t>o is defined. For any Banach space 

H, denote by L 2 r (0,T;ff) the Banach space consisting of all if -valued {J^jj^o-adapted processes 

X(-) such that E(|X(-)| 2 2 , 0T . H J < oo, by Lj?(0,T;H) the Banach space consisting of all if -valued 



{■^t}t>o-adapted bounded processes, and by Lp(£l;C([0,T];H)) the Banach space consisting of 



all if -valued {J^ji^o-adapted continuous processes X(-) such that E(|X(-)| 2 7 ,r 0T , H ^) < oo. All of 

these spaces are endowed with the canonical norm (Similarly, one can define Ljr(Q;C([0,T];H)) 
for any positive integer k). 

Throughout this paper, we assume that the functions b 1 ^ G C 1 (G) (i,j = 1,2, ••■ , n) satisfy 
5*7 = \y> 1 a nd, for some constant sq > 0, 

n 

J2 b ij ee > soiei 2 , v(x,o e g x R n . (1.1) 



*This work is partially supported by the NSF of China under grants 10831007 and 60974035 and the Grant 
MTM2011-29306 of the MICINN, Spain. 

'School of Mathematical Sciences, University of Electronic Science and Technology of China, Chengdu 610054, 
China; and BCAM-Basque Center for Applied Mathematics, Alameda de Mazarredo, 14, E-48009 Bilbao, Basque 
Country, Spain. E-mail: luqi59@163.com. 

'School of Mathematics, Sichuan University, Chengdu 610064, China. E-mail: zhangjcu@scu.edu.cn. 



Consider the following stochastic hyperbolic equation: 



dz t 



n 



{b ij z Xi ) x At = (hzt + b 2 -Vz + b 3 z) dt + (b 4 z + g)dB(t) in Q, 







z(0) = zq, z t (0) = z 1 



on E, 
inG. 



;i.2) 



Here, zj = ^|, z Zi = ^-, and bi (1 < i < A) are some suitable known functions to be given later; 
while {zq,z\) G L 2 (n,F ,P;H^(G) x L 2 {G)) and 5 G L 2 T (0,T;L 2 (G)) are unknown. Physically, 5 
stands for the intensity of a random force of the white noise type. Put 



H T ^L 2 T (n;C([Q,T];H l (G)))f^L 2 An;C\[Q,T];L 2 (G))). 



;i.s) 



It is clear that Ht is a Banach space with the canonical norm. Under suitable assumptions (the 
assumptions in this paper are enough), for any given (zo, z±) and g, one can show that the equation 
(1.2) admits one and only one solution z = z(zo,zi,g)(t,x,u>) G Ht (see [12]). We will also denote 
by z(zo,zi,g) or z(zo, z±,g)(t) the solution of (1.2). 

In this article, the random force L gdB is assumed to cause the random vibration starting 
from some initial state (zq,zi). Roughly speaking, our aim is to determine the unknown random 
force intensity g and the unknown initial displacement zq and initial velocity z\ from the (partial) 
boundary observation ^L nT \ r and the measurement on the terminal displacement z(T), where 
v = v{x) denotes the unit outer normal vector of G at x G T, and Tq is a suitable open subset 
(to be specified later) of T. More precisely, we are concerned with the following global uniqueness 
problem: Do ]^;(zo, z ii9)\(nT)xT = ^ an< ^ z ( z 0' z i'5')(^) = in G, P-a.s. imply that g = in Q 
and zq = z\ = in G, P-a.s. ? 

In the deterministic setting, there exist numerous literatures addressing the inverse problem of 
PDEs (See [5, 8] and the rich references cited therein). A typical deterministic inverse problem 
close to the above one is as follows: Fix suitable known functions a(-,-) and /i(-,-) satisfying 
min |/i(t,x)| > 0, and consider the following hyperbolic equation: 

(t,x)€Q 



ztt-Az = a(t,x)z + fi(t,x)f 2 (x) in Q, 

z = on S, 

k z(0) = 0, z t (0) = z-y in G. 



(1.4) 



In (1.4), both Z\ and ji are unknown and one expects to determine them through the boundary 
observation ^ c Tl r ■ As shown in [11], by assuming suitable regularity on functions a(-, •),/«(-, •) 
(i = 1,2) and zi(-), and using the following transformation 



, . d ( z(t,x) 
y = y(t,x) = — 



(1.5) 



dt \fi(t,x) 
this inverse problem can be reduced to deriving the so-called observability for the following wave 



equation with memory 

y tt - Ay = oij/t + a 2 • Vy + a 3 y 



+ / [ci(t, s, x)y(s,x) + c 2 (t, s,x) ■ Vy(s,x)]ds in Q, 
Jo 



y = 

y(o,x) = ^M., w ( 0> x) = / 2 (x)-^(W ZlW 



A(0,x) 



|/l(0,x)P 



on £, 
inG, 



where aj(-, •) (i = 1, 2, 3) and c$(-, -, •) (i = 1, 2) are suitable functions. Concerning this problem, if 
z\ is known and both functions a(-, ■) and /i(-, •) are independent of the space variable x, i.e., there 
is only one unknown in (1.4), then the corresponding inverse problem is now well-understood (e.g. 
[9, 10] and the references therein). The main tool in the later case is to use the Duhamel principle, 
instead of the transform (1.5), to reduce the problem to the observability estimate for some wave 
equation (without memory). 

Stochastic partial differential equations (PDEs for short) are used to describe a lot of random 
phenomena appeared in physics, chemistry, biology, control theory and so on. In many situations, 
stochastic PDEs are more realistic mathematical models than the deterministic ones. Nevertheless, 
compared to the deterministic setting, there exist a very limited works addressing inverse problems 
for stochastic PDEs. In this respect, we mention [1] for a study of an inverse medium scattering 
problem for the random Helmholtz equation. Also, we refer to [2, 4] for several results on the 
estimation problems of some random and stochastic PDEs when the noise intensity tends to zero. 
To the best of our knowledge, there is no paper considering the inverse problem for stochastic 
hyperbolic equations. 

One may meet substantially new difficulties in the study of some inverse problems for stochastic 
PDEs. For instance, unlike the deterministic PDEs, the solution of a stochastic PDE is usually 
non-differentiable with respect to the variable with noise (say, the time variable considered in 
this paper). Also, the usual compactness embedding result does not remain true for the solution 
spaces related to stochastic PDEs. These new phenomenons lead that some effective methods for 
solving inverse problems for deterministic PDEs (see [9] for example) cannot be used to solve the 
corresponding inverse problems in the stochastic setting. Especially, one can see that none of the 
methods for solving the above inverse problem for the equation (1.4) can be easily adopted to solve 
our inverse problem for the stochastic hyperbolic equation (1.2), even if g is assumed to be of the 
form 

g(t,x,uj) = gi(t,w)g 2 (x), V (t,x,u) G Q x fl, (1.6) 

with a known nonzero stochastic process <7i(-, ■) £ Ljr(0,T) and an unknown deterministic function 
<72(") £ L 2 (G). For these reasons, it is necessary to develop new methodology and technique for 
treating inverse problems for stochastic PDEs. 

In this paper, we will use a global Carleman estimate to solve the above formulated inverse 
problem for the equation (1.2). As far as we know, [12] is the only published reference addressing 
the Carleman estimate for stochastic hyperbolic equations. In [12], under suitable assumptions, 
the following estimate was proved for the solution z of (1.2): 



\(z(T),z t (T))\ L 2( a:FT! p. H i( G ) xL 2( G ^ 



<C 



dz 



dv 



L%(0,T;L2(r )) 



+ \9\l%(0,T;L 2 (G)) 



;i.7) 



(Here and henceforth, C is a generic positive constant, depending only on T, G, Tq and sq, which 
may be different from one place to another). Noting however that the (random) source g appears 
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in the right hand side of (1.7), and therefore, the estimate obtained in [12] does not apply to the 
inverse problem considered in this work. In order to solve our stochastic inverse problem, we have 
to establish a new Carleman estimate for (1.2) so that the source term g can be bounded above 
by the observed data. Hence, we need to avoid employing the usual energy estimate because, when 
applying this sort of estimate to (1.2), the source term g would appear as a bad term. Meanwhile, 
since we are also expected to identity the initial data, we need to bound above the initial data by the 
observed data, too. Because of this, we need to obtain the estimate on the initial data and source 
term in the Carleman inequality simultaneously. Therefore we cannot use the usual "Carleman 
estimate" + "energy estimate" method (which works well for the deterministic wave equation, see 
[3]) to derive the desired estimates. This is the main difficulty that we need to overcome in this 
paper. 

The rest of this paper is organized as follows. In Section 2, we state the main result of this paper. 
Some preliminary results are collected in Section 3. Finally, Section 4 is addressed to proving the 
main result. 

2 Statement of the main result 

To begin with, we introduce the following conditions: 

Condition 2.1 There exists a positive function d(-) £ C 2 (G) satisfying the following: 

1) For some constant no > 0, it holds that 

n n n 

E{ E [^ ij '(b^d Xil ) Xjf -b^'d Xil ]}ee >m<,i>w, 



i,j=l i',j'=l i,j=l 



(2.1) 



V(*,£V--,neGx 



2) There is no critical point of d(-) in G, i.e., 



min|Vd(x)| > 0. (2.2) 

xeG 

Remark 2.1 If (6' J ) i<ij< n is the identity matrix, then, by taking d{x) = \x — xo| 2 withxQ ^ G, one 
sees that Condition 2.1 is satisfied. Condition 2.1 was introduced in [3] to show the observability 
estimate for hyperbolic equations. We refer to [3] for more explanation on Condition 2.1 and 
illustrative examples. 

It is easy to see that if d(-) G C 2 (G) satisfies Condition 2.1, then for any given constants a > 1 
and dsl, the function d = ad + b still satisfies Condition 2.1 with jjlq replaced by a/io- Therefore 
we may choose no as large as we need in Condition 2.1. Now we choose < cq < c± < 1, no > 4 
and T satisfying the following condition: 

Condition 2.2 

1) ^o - 4ci - c > 0, 



n 



2 ) (8c +c) £ biJdxid *i > 4C?T2 > £ h%HxA 

i,j=l i,j=l 



Remark 2.2 Since N^ b l ^d Xi d Xj > 0, it is easy to see that one can always choose hq in Condi- 

tion 2.1 large enough so that Condition 2.2 holds true. We put it here simply to emphasize the 
relationship among < cq < C\ < 1, /Zq > 4 and T . 



In the sequel, we choose 



To = <, •- 



GT: J2b^d Xi ( X y(x)>0]. 



(2.3) 



Also, we assume that 

6i G L^(Q,T;L°°(G)), b 2 G L^(0,T; L°°(G;M n )), 

b 3 G L^(0,T;L n (G)), 64 G L^(0,T; L°°(G)). 

In what follows, we use the notation: 



^A„ |2 



l&l 



L??(0,T;L°°(G)) 



+ 1^2 



+ 16; 



2lLJ?(0,T;L°°(G;R™)) "•" \ u 3\L^(0,T;L n (G)) "•" r4|L£?(o,T;L°°(G)) 



+ I&4 



+ 1. 



(2.4) 



(2.5) 



The main result of this paper can be stated as follows. 



Theorem 2.1 Let b{ (1 < i < 4) satisfy (2.4), and let Conditions 2.1 and 2.2 hold. Assume that 
the solution z G Ht of (1.2) satisfies that & = on (0, T) x To and z(T) = in G, P-a.s. Then 
g = in Q and Zq = Z\ = in G, P-a.s. 

Several remarks are in order. 

Remark 2.3 Similar to the inverse problem for (1-4), and stimulated by Theorem 2.1, it seems 
natural and reasonable to expect a similar uniqueness result for the following equation 



dzt - Yl (b ij z Xi ) Xj dt = {bizt + 6 2 • Vz + b 3 z + f)dt + b A zdB(t) in Q, 



z = 
k 2(0) = zq, zt(Q) = zi 



on S, 
in G, 



(2.6) 



in which zq, z\ and f are unknown and one expects to determine them through the boundary 
observation ^\( 0T \ xT an< ^ ^ e ^ er ™ na ^ measurement z(T). However the same conclusion as that 
in Theorem 2.1 does NOT hold true even for the deterministic wave equation. Indeed, we choose any 
y G Cq°(Q) so that it does not vanish in some proper nonempty subdomain of Q. Put f = utt — An. 
Then, it is easy to see that y solves the following wave equation 

y tt -Ay = f in Q, 

y = 0, on S, 

I 1/(0) = 0, y t (0) = mG. 

One can show that y(T) = in G and ^ = on S. However, it is clear that f does not vanish 
in Q. This counterexample shows that the formulation of the stochastic inverse problem may differ 
considerably from its deterministic counterpart. 



Remark 2.4 From the computational point of view, it is quite interesting to study the following 
stability problem (for the inverse stochastic hyperbolic equation (1.2)): Is the map 



dz 

-^-(z ,zi,g) 

dv 



(o,T)xr 



x z(zQ,z 1 ,g)(T) — > (zo,zi,g) 



continuous in some suitable Hilbert spaces? Unfortunately, we are not able to prove this stability 
result at this moment. Instead, from the proof of Theorem 2.1 (See Theorem 4-1 in Section 4), it 
is easy to show the following partial stability result, i.e., for any solution z £ Ht of the equation 
(1.2) satisfying z(T) = in G, P-a.s., it holds that 



dz 



dv 



L2.(0,T;L 2 (r )) 



\{ z Oi z l)\L 2 (n,To,P;Hi,(G)xL 2 (G)) + \vT^t9\L%(0,T;L 2 (G)) ^ C 

Especially, if g is of the form (1.6) (with <7i(-,-) £ Ljr(0, T) \ {0} and <72(") £ L 2 (G)), then the 

dz 



following estimate holds 

\{ z 0-, z l)\L 2 (n,T ,P;Hl(G)xL 2 {G)) + Ml 2 (G) < C 



dl 



L%(0,T;L*(r )) 



Remark 2.5 The inverse problem considered in this work is a sort of inverse source problems. It 
would be quite interesting to study the global uniqueness and stability of inverse coefficient problems 
for stochastic PDEs but this remains to be done, and it seems to be a very difficult problem. 

Remark 2.6 It is also interesting to study the same inverse problems but for other stochastic 
PDEs, say the stochastic parabolic equation, the stochastic Schodinger equation, the stochastic plate 
equation and so on. However, it seems that the technique developed in this paper cannot be applied 
to these equations. 



3 Some preliminaries 

In this section, we collect some preliminaries which will be used later. 

First, we show the following hidden regularity result for the solution z to the equation (1.2) (This 



result means that the observation of the normal derivative of z makes sense, i.e., \^;\l 2 (o,t ; l 2 (T )) 
+oo). 



< 



Proposition 3.1 Let b{ (1 < i < 4) satisfy (2.4)- Then, for any solution of the equation (1.2), it 
holds that 



dz 



dv 



L 2 (0,T;L 2 (T)) 



< e 



CA 



K z 0,^l)lL2(Q i jr Oi p. H i( G ) xL 2( G )) + Ml2.(o,T;L 2 (G)) 



(3.1) 



Remark 3.1 In [12], the author proved Proposition 3.1 when (6 J - , )i<i ) j<n is an identity matrix. 
The proof of Proposition 3.1 for the general coefficient matrix (b l:) )i<ij<n is similar, and therefore 
we give below only a sketch of the proof. 



Proof of Proposition 3.1: Since V G C 2 , one can find a vector field h = (/i 1 ,--- ,h n ) £ 
C 1 (W n ;W n ) such that h = v on T (see [6]). A direct computation shows that 

n n n 



J2 [2(h ■ Vz) Y, b l3 z Xj + h* (z 2 - £ b ij z Xi ; 



dt 



i=i 



i=i 



*j=i 



n n 

2 { [ dzt ~ E (^'^i)^] fc ■ Vz - d ^ h ■ Vz ) - E b ij z Xi z Xk h k x . } 



dt. 



(3.2) 



*,i=i 



i,j,k=l 



—z t divhdt + \ z x z Xi dw (b tJ h)dt. 

Integrating the identity (3.2) in Q, taking expectation in O and using integration by parts, we 
obtain the inequality (3.1) immediately. □ 

Next, we recall the following known result. 

Lemma 3.1 ([12]) Let p v G (^((O.T) x R n ) satisfy 

p^=p>\ i,j = 1,2,- ■■ ,n, (3.3) 

£, f, ty € C 2 ((0,T) x M n ). Assume that u is an H 2 oc (M. n ) -valued {J- t}t>o~ adapted process such that 
ut is an L 2 (M. n ) -valued semimartingale. Set 6 = e e and v = 6u. Then, for a.e. a; € W 1 and P-a.s. 

n n 

e(- 2£ t v t + 2 J2 P iJ ^v Xj + *v) \du t - Y, (P i] u Xt ) Xj dt 



*j=i 



»j=i 



n n 



+ E [ E {^P l ' f ^v^ Xj/ -p^p^'l Xi v Xif v Xj ,)-2fH t v Xi v t +p'H Xi v 2 



i,j=l i',j'=l 



+Vp ij v Xi v- [A£ Xi + 



* 



' W 



dt 



*/ 



+d[ Y P iJ£ tv Xi v Xj - 2 Y J j ZxiV Xj vt + £ t v 2 - *v t v + [A£ t + ^-)r 2 

i,j=l i,j=l 



(3.4) 



n n 

{ [in + E (P^'k " *] v * " 2 E [&%)t +P ij ttj]v x , 



»J=1 



1'f 



«J=1 



n n 



+ \ I fn J J 



fxiVx,- 



z',j'=l 



+Bv 2 + ( - 2f m + 2 £ p%^. + #v) }dt + 9 2 £ t (du t ) 2 , 



»,j=i 



where (dut) 2 denotes the quadratic variation process of ' ut, A and B are stated as follows: 



A = (£ 2 - £ tt ) - J2 (P iJ ^£ Xj ~ P%^ ~ P ij i XiXj ) ~ % 

n 1 n 

B^AV + (A£ t ) t - Y {Ap ij £ Xi ) Xj + x [*» - £ (p^)*. 



(3.5) 



»,j=i 



*j=i 



4 Proof of the main result 

This section is devoted to proving Theorem 2.1. As mentioned before, we will prove Theorem 2.1 
by establishing a new Carleman estimate for the equation (1.2). 
In the sequel, we choose 

9 = e e , e = \[d(x)-a(t-T) 2 ], 

where A > is a parameter, d(-) is the function given in Condition 2.1, and c\ is the constant in 
Condition 2.2. 

Our global Carleman estimate for (1.2) is as follows. 

Theorem 4.1 Let hi (1 < i < 4) satisfy (2.4), and let Conditions 2.1 and 2.2 hold. Then, there 
exists a constant A > such that for any A > A and any solution z £ H? of the equation (1.2) 
satisfying z{T) = in G, P-a.s., it holds that 



E I 2 (A|*i| 2 + A|Vz | 2 + A 3 |z | 2 )^ + AE I (T - t)9 2 g 2 dxdt 
Jg Jo 



< CXE f 

Jo J To 



dz 



2 

dTdt. 



(4.1) 



dv 

Proof: In what follows, we shall apply Lemma 3.1 to the equation (1.2) with 

n 

u = z, p«=&tf, y = itt+Y,(b i n Xi ) Xj -\c , 

(Recall Condition 2.2 for the constant Co), and then estimate the terms in (3.4) one by one. 

In the sequel, for A > 0, we use 0(X r ) to denote a function of order A r for large A. The proof 
is divided into three steps. 

Step 1. In this step, we analyze the terms which stand for the "energy" of the solution of (1.2). 
As the Carleman estimate for deterministic partial differential equation, the point is to compute 
the order of A in the coefficients of \vt\ 2 , \X7v\ 2 and \v\ 2 . Since the computation is very close to that 
in the proof of [7, Theorem 1.2.1], we give here only a sketch. 

First, it is clear that the coefficient of \vt\ 2 reads: 



n 

itt ' ' 



+ E ( 6 ^*i)*i - * = Ac o- ( 4 - 2 ) 



Further, noting that b %3 (1 < i,j < n) are independent of t and £t Xj = %-xjt = 0) we fi n d that 

n 

Y,( biii *i)t + b ij £ tXj v Xi v t = 0. (4.3) 

Further, by Condition 2.1, we see that 

n n 

E {(^)t+ J2 [^'{U'n^) X] ,-{b^h^'£ Xi ,) X] ] +Wi}v Xi v.. 



i,j=l i',j'=l 






(4.4) 
> A(/x - 4ci - c ) E b lj v Xi Va 



"3 



Further, in order to compute the coefficient B of \v\ 2 , recalling (3.5), we find that 



A = X 2 



4c 2 (t _ T) 2_ J2 b^d Xi d Xj ] + O(A). 



Hence, by the definition of B (in (3.5)), we conclude that 

n n n 

B = (4 Cl + c ) E ^^M^ 3 + E Y, bt3 d Xl (.V' j 'd Xi ,d Xj ,) Xj \ 3 

-4(8c? + c cl)(t - T) 2 A 3 + 0(A 2 ). 
Recall the following estimate in [7]: 



n n 



W^ ^d Xi d Xj < E E bt3 d Xl (b t,f d Xif d Xjl ) 
i,j=l i,j=li',j'=l 



Therefore, by Condition 2.2, we obtain that 



(4.5) 



(4.6) 



(4.7) 



B > (4ci + c ) E b ij d Xi d Xj X 3 + M o E bl3 d Xl d Xj \ 3 

i,j=l i,j=l 

-4(8ci + co)cf(t - T) 2 A 3 + 0(A 2 ) 
= (4ci + c ) E b ij d Xi d Xj \ 3 + 0(A 2 ). 

Hence, there exists a Ao > such that for any A > Ao, it holds that 

5t> 2 > CAV. 



(4i 



(4.9) 



Step 2. In this step, we analyze the terms corresponding to t = and t = T. For the time 
£ = 0, we have 



n n 1 

E b ij £ t v Xt v Xj - 2 E b ij £ Xi v Xj v t + £ t v\ - ^v t v + Ul t + -^ 



2 Cl TA E b ij v Xi v Xj - 2A E ^ «*,-«* - a( - 2 Cl + E (^4J*. 



»j=i 



*,i=i 



C f if 



*J=1 



2ciTAw t 2 + [2ciT(4c 2 T 2 - ^ b ij d Xt d Xj )\ 3 + 0(A 



»j=i 



> 2 Cl TA J^ 6^.^. - A( E ^4,4,) 2 E ^;^ - A( E ^'<^.J V 

i,j=i *,i=i *,i=i *,i=i 



i n 



+2ciTAw t 2 - t> t 2 + 



2 Cl T(4c 2 T 2 - J2 b ij d Xi d X] j\ 3 + 0(A 2 



(4.10) 



By Condition (2.2), it follows that 



4c 2 T 2 - Y b^d Xi d Xj > 
and that 

n i 

2ciT-^b ij d Xi d Xj y >o. 

Hence there exists a Ai > such that for any A > Ai and when t = 0, it holds that 

n n -. 

Y b ij £tv Xi v Xj - 2 Y b ij £ Xi v Xj v t + £ t v 2 - Vv t v + (A£ t + -* t 



> C 



i,j=l 



+ |Vv| 2 ) + AV|. 



On the other hand, since £t{T) = 0, for t = T, it holds that 

n ra 

Y biJ ttv Xi v Xj - 2 J^ b ij £ Xi v Xj v t + £ t v 2 - Vv t v + (Alt + ~«, 



(4.11) 



t v- 



*j=i 



»j=i 



(4.12) 



= -2 ^ b ij £ Xi v Xj v t - Vv t v. 

Noting that z(T) = in G, P-a.s., we have f(T) = and v x .(T) = in G (j = 1,2,-- ,n), P-a.s. 
Thus, from the equality (4.12), we end up with 



{n n 1 I 

Y b ij £ t v Xi v Xj - 2 Y b ij £ Xl v X] v t + £ t v 2 - Vv t v + (Al t + -* t )v 2 \ 



0, P-a.s. (4.13) 



t=T 



Step 3. Integrating (3.4) in Q, taking expectation in f2 and by the argument above, for 
A > max{Ao, Ai}, we obtain that 



„ n n 

E / e{ ( - 2€m + 2 ^ ^4^ + *«) [dz t - Y (b ij z Xi ) Xj dt\ ]dx 

„ n n 

+ AE / E E (2^& iY 4 i ,« Xi t; a;j ,-6«6^''d a;i ^. / t; Xj ,) I / i ,-dS 



(4.14) 



>C^E 



Q 



Az 2 + A|Vz| 2 + A 3 * 2 J + ( - 2^ + 2 J^ VU x .v Xi + *7 



dxdt 



+E 



(V 



3|„ |2 



A(|vz r + Nir) + A d |z 



dx + K I 6% 
Q 



(dz t ) 2 }. 



For the boundary term, noting that z = on S, it is easy to show that 

„ n n 

E / E E (w i V' i '<^ i .Vx i v Xj ,-V i V' i 'd Xi v x . l v Xj ,y j d2 

„ n n 

= e (Y bii wi){ E hi ' j ' d ^ 



(4.15) 
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*J=1 

dv 2 
dv 



From inequality (4.14) and equality (4.15), we obtain that 

,. n n 

E / e[ ( - 2^ t + 2 J^ ^4^- + *«) [d«t - J^ (6^)*,^] }cte 

„ n n 

+AE / ( J>^/)( J] ^"^ 

n 

^Az 2 + A|Vz| 2 + A 3 z 2 ) + ( - 2£ t v t + 2 J^ ^4^- + * 

+E / 2 [a(|Vzo| 2 + |^i| 2 ) + A 3 |z o | 2 1^ + AE [ (T-t)9 2 (b 4 z + g) 2 dxdt} 
Jg l j 7q j 

By means of 



>C^E 



dxdt 



b A z + gf > -g 2 - 2b\z\ 



we get 



AE / (T - £)6> 2 (& 4 ,z + #) 2 d:rcft > -AE / (T - t)6 2 g 2 dxdt - 2ATE / Q 2 b\z 2 dxdt. 
JQ 2 JQ JQ 



'Q z JQ 

On the other hand, by equation (1.2), it is clear that 



E f 0[ ( - 2^ + 2 £ ^4^ + *«) [dzt - £ (b lj z x J Xj dt\ }dx 

< E f ( - 2^ + £ 6°'4i^ + Vvfdxdt + c||6i| 2 ?(0T;LOO(G)) E / 2 z 2 dxdt 

^2|l^(o,T;L°°(G,IR")) + P3|l^(o,T;L"(G)) 



E / # 2 |Vz| 2 da;(ii 
Q 



+ ^ 2 |Ml°°(o,T;L"(G)) E ' ® 2z2 dxdl 



(4.16) 



(4.17) 



(4.18) 



Finally, taking A = max {CA, A , Ai}, combining (2.3), (4.16), (4.17) and (4.18), for any A > A, 
we conclude the desired estimate (4.1). □ 

We are now in a position to prove Theorem 2.1. 

Proof of Theorem 2.1 : Since ^ = on So, P-a.s., we know the right hand side of inequality 
(4.1) is zero. Therefore, it follows that 



and that 



E / 9\\\zi\ A + A|Vz r + y\z \ z )dx = 
'g 



E / (T - t)6 2 g 2 dxdt = 0. 
JQ 



(4.19) 



(4.20) 



From the equality (4.19), we find zq = z\ = in G, P-a.s. By means of the equality (4.20), we see 
g = in Q, P-a.s. □ 
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